Entanglement spectrum of one-dimensional extended Bose-Hubbard models 



O 
(N 



C/2 



o 

(N 
> 

in 

^' 

o 



X 



Xiaolong Deng and Luis Santos 
Institut fiir Theoretische Physik, Leibniz Universitdt Hannover, Appelstr. 2, 30167 Hannover, Germany 

(Dated: January 20, 2013) 

The entanglement spectrum provides crucial information about correlated quantum systems. We 
show that the study of the block-like nature of the reduced density matrix in number sectors and the 
partition dependence of the spectrum in finite systems leads to interesting unexpected insights, which 
we illustrate for the case of a ID extended Hubbard model. We show that block symmetry provides 
an intuitive understanding of the spectral double degeneracy of the Haldane-insulator, which is 
remarkably maintained at low on-site interaction, where triple or higher site occupation is significant 
and particle-hole symmetry is broken. Moreover, surprisingly, the partition dependence of the 
spectral degeneracy in the Haldane-isulator, and of a partial degeneracy in the Mott-insulator, are 
directly linked to the, in principle unrelated, density-density correlations, and presents an intriguing 
periodic behavior in superfluid and supersolid phases. 



Ultra-cold gases in optical lattices [H [3 allow for the 
controlled study of Hubbard- like models Q , fundamen- 
tal for the description of strongly correlated systems. Al- 
though until recently only on-site interactions have been 
relevant, inter-site interactions may play a crucial role in 
new experiments on gases with dipolar interactions [bi- 
polar lattice gases may allow for the simulation of ex- 
tended Hubbard models, which play an important role in 
e.g. the analysis of high- Tc superconductivity Q, being 
the simplest models that capture the interplay between 
strong correlations and charge-ordering effects. 

Most remarkably, it has been recently shown that ID 
polar lattice bosons, in addition to conventional locally- 
ordered phases, may allow for the so-called Haldane- 
insulator (HI) phase (a Haldane-liquid phase has 
been recently proposed for polar fermions Q). The HI 
lacks local order, but presents non-local string-like or- 
der [13], resembling to a large extent the Haldane phase 
in integer spin chains [llj. The Haldane phase is one of 
the paradigms of topological quantum phases, which due 
to the lack of local order parameter fall beyond the Lan- 
dau paradigm of condensed-matter physics. Topological 
phases have attracted a growing interest in recent years, 
most relevantly their entanglement properties. In par- 
ticular, topological phases present an anomalous entan- 
glement entropy, characterized by a topology-dependent 



^ ] universal additive constant [12, llSj 



The entanglement entropy is defined as the von Neu- 
mann entropy, S — —TrpA^ogpA, associated to the re- 
duced density matrix, pA , obtained after partitioning the 
system into two parts, A and B, and tracing out B. How- 
ever, Pa contains much more information than its mere 
entanglement entropy. In particular, it has been shown 
that the set of eigenvalues of pA, the so-called entan- 
glement spectrum (ES), provides key insights concerning 
quantum correlations in many-body systems (l4| . The 
intriguing properties of the ES in different physical sys - 
tems have attracted recently a major interest [15l - ll8j . 
Interestingly, it has been shown that the ES of the Hal- 
dane phase in ID integer-spin chains is characterized by a 



double degeneracy protected by a set of symmetries [19| . 
Moreover, the difference between the two largest eigen- 
values of the ES (Schmidt gap) presents at the boundary 
of the Haldane phase in ID 5 = 1 spin chains scaling 
properties which coincide to those of the corres pon ding 
local order parameter in non-topological phases [20j. 

In this Letter, we consider the ID extended Bose- 
Hubbard model (EBHM) to illustrate the interesting, 
and somewhat unexpected, insights provided by a careful 
study of the ES. As for the Haldane phase, the HI phase 
is characterized by a doubly-degenerated ES. Our numer- 
ical results show that, remarkably, this double degener- 
acy (as suggested in i^M,) is maintained at low on-site re- 
pulsions, where significant triple, or higher, site occupa- 
tion breaks particle- hole symmetry. We show that, as re- 
cently discussed [3,[2l[, the analysis of the block structure 
of PA may offer interesting insights. In particular, the 
block symmetry of pA permits an intuitive understand- 
ing of the doubly-degenerated ES of the HI, unveiling 
partial double degeneracy for the Mott-insulator (MI). 
Moreover, we show that the partition-dependence of the 
ES in finite systems, which to the best of our knowl- 
edge was not yet considered, provides interesting sur- 
prising insights on the relation between ES and excita- 
tion spectrum. In particular, the Schmidt gap in the 
HI, and a similar one for the MI, present an exponential 
decay in its partition-dependence, with a characteristic 
length which, remarkably, is identical to that of the, a 
priori unrelated, density-density correlations. Moreover, 
the partition-dependence provides interesting informa- 
tion about other phases, presenting an intriguing peri- 
odicity in superfluid (SF) and supersolid (SS) phases. 

We consider below the ID EBHM 
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where &|, bi are creation and annihiliation operators for 
bosons at site i, rii — bjbi, t is the hopping amplitude. 



2 



and U and V are, respectively, the coupling constants for 
on-site and inter-site interactions. As mentioned above, 
inter-site interactions appear naturally for polar particles 
in lattices. Although for polar lattice gases interactions 
beyond nearest neighbors may play a role at low energies, 
model (dJ contains already most of the relevant physics. 

The Holstein-Primakoff transformation — {2S — 
n,f'%,, S- = bl{2S - mfl'^ and Sf = S - n, links the 
EBHM with the anisotropic Heisenberg spin-S model: 
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where J^/J = VS/2t and D/J = US/At characterize, 
respectively, the exchange and single-ion anisotropy. As- 
suming (5"^) — (i.e. (rij) = S), both models are equiva- 
lent if 5ni ^ Hi — S S and the maximal site occupation 
nmax < 25'. Below we consider (rii) = 1 which relates to 
S — 1. However, since S = 1 does not fulfill the afore- 
mentioned criteria significant departures between both 
models may be expected, especially for low U/t. 

For D/J > and JJJ > 0, model ^ presents 
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three possible phases 
Sz — is favoured 
large-D phase, where the 
limi,-_,i^_(-l)l*-J' 
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For sufficiently large D, 
and the system enters the 
Neel order O^'f^^^''' = 
{S^S") = 0, and the non-local string 
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tween the Neel and large-D phases, the Haldane phase 
is found, a gapped phase with antiferromagnetic spin or- 
der (Oaring 7^ 0) but no Spatial order (OnocI — 0)- 

At large-enough U/t, the EBHM phases for {m) = 

1 map those of = 1 chains. The equivalent of 
the large-D phase is a MI with one boson per site, 
with string order parameter 6] Ogtring = liiiir-i.oo((l — 
no)e*'^5i;fc=o"fc(l — n^)) — 0, density-density order param- 
eter Cnn = liiiir-ioo[('T'0'^r) — ('^o)("-r)] = 0, and parity 
order parameter [7| Oparity = limr_>.oo(e"^^'^=o"'=) ^ 0. 
A density wave (CDW) of alternate sites with and 

2 bosons is the counterpart of the Neel phase, with 
Ostring 7^ 0, Oparity ^ 0, and Cnn 0. The equivalent 
of the Haldane phase, the HI [1], has Oparity = Cnn = 0, 
Ostring 7^ 0. For low U/t, the system may be SF or SS, 
with polynomically decaying single particle correlation. 

We investigate the EBHM using density-matrix renor- 
malization group (DMRG) [22| with up to L = 100 sites, 
keeping ~ 300 optimal states. Since the EBHM is num- 
ber conserving, we implement our DMRG for sectors 
with a specific boson number. We consider maximally 4 
bosons per site, which we have checked to be enough for 
U/t > 0.5. A proper investigation of the bulk properties 
of the HI using open boundary conditions requires avoid- 
ing the spurious effects introduced by edge states p3| . 



We have employed two independent ways of handling 
edge states, ensuring that our results are method inde- 
pendent. A first method consists in placing N = L + 1 
bosons, resulting in a magnetization Sf = —1 which, 
due to the Haldane gap, is accomodated at the edges, 
5f = 5£ = —1/2, with no energy cost. The edge states 
are then polarized. An alternative method consists in 
adding a site with a hard-core boson at each chain bound- 
ary (sites and L + 1). This method is equivalent to that 
employed in S* = 1 in Ref. By properly tuning the 
coupling between the bulk and the extra sites, we may 
achieve (no) = (rii+i) sa 1/2, i.e. the same probability 
of effective spin-1/2 up and down states. The spin-1/2 
forms a singlet with the edge spins, destroying the edge 
states. Note, however, that the incommesurability as- 
sociated with the handling of the edge states leads to 
a commesurate-incommesurate transition of the MI into 
SF. The analysis of the MI phase demands hence L — N. 

We recall that the ES is obtained after dividing the 
system into two parts, A (with La sites) and B, trac- 
ing out B, and diagonalizing the reduced density matrix 
Pa- Note that due to number conservation in the whole 
system, pA is block diagonal, with blocks p{Sa) char- 
acterized by an effective magnetization Sa = La — Na 
given by the number of particles Na in A. The ES is 
hence grouped in sectors with a given S^, Xi{S^), with i 
the index of the eigenvalue (in decreasing order). The 
structure of Ai(5^) (Figs. [IJ offers an interesting in- 
sight into the degeneracies of the ES. Fig. [lla) shows 
that, as for the Haldane phase of the S — 1 chain (l9i |. 
a doubly degenerate ES is recovered for the HI phase. 
Fig. [IJa) shows that the double degeneracy is directly 
linked to the symmetry p{Sj^ + s) = p(S\ — s — 1), with 
s = 0, 1, . . . [25], which may be linked to number conser- 
vation and bond-centered symmetry fl9| . Remarkably, 
we observe a basically perfect double degeneracy in the 
HI phases even at U/t as low as 1, where occupations 
rii > 2 are very significant, hence clearly breaking parti- 
cle/hole symmetry. Finally note that {Ai(5^)} present a 
Gaussian dependence around S^- = —1/2, resembling the 
ES of permutationally-symmetric models jl^l, and that 
blocks with 5'|_^> 1(< —2) appear due to solitons in the 



string-order [26[ and sites with Ui > 2. 



Interestingly, the ES of the MI (Fig. [1Kb)) presents 
block symmetry as well, p{Sa + s) = p(Sa — s). How- 
ever, the eigenvalues of S*^ = are not degenerated. As 
a result, the spectrum of the MI presents a mixture be- 
tween a partial double degeneracy and non-degenerate 
eigenvalues. In particular, a gap is observed between 
the two largest eigenvalues of the spectrum (Schmidt 
gap [2^), whereas the second and third eigenvalues are 
degenerated. Finally, in the CDW the block symmetry is 
lost (Fig. [1] (c)), and hence the ES does not present any 
systematic degeneracy. All cases show an approximate 
Gaussian dependence of , which may be employed 
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FIG. 1: Largest eigenvalues of the ES for different S% 
sectors, for U/t = 5, L = 100, and V/t = 3.3 (HI) (a), 2.0 (MI) 
(b) and 4.0 (CDW) (c). The results for X^^ are compared 
with a Gaussian distribution (dotted). 



for designing more efficient DMRG codes. 

In finite systems, the ES depends on the position of 
the partition, i.e. on La- Interestingly, as we show 
in the following, this La dependence provides interest- 
ing non-trivial information about both locally and non- 
locally ordered phases. We plot in Fig. ^a) the Schmidt 
gap A{La) = Ai — A2 between the two largest eigen- 
values of Pa for different values of La, starting from 
La — 1. Sufficiently away from the chain borders, we 
find that in the HI A{La) and also all X2j-i — X2j de- 
cay exponentially. Remarkably, as shown in Fig. [Hb), 
the typical length characterizing this exponential decay, 
^E, is identical to the correlation length, ^„„, associ- 
ated to the exponentially decaying density-density corre- 
lation in the bulk C„„(i) = {njUj+i) — {nj){nj+i). Since 
^nn is inversely proportional to the Haldane-like gap for 
number-conserving excitations (inset in Fig. [ll^b)), our 
results show that there is a non-trivial, and somewhat 
unexpected, link between the partition dependence of 
the ES and the energy gap characterizing the excita- 
tion spectrum. Moreover, a similar analysis for the La 
dependence of the ES may be performed for the gap 
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FIG. 2: (a) Partition dependence of the Schmidt gap (and 
As - As) for L = 150, U/t = 5, and V/t = 3.2 (HI), com- 
pared to the density-density correlation C„„(i); (b) Decay 
length of the Schmidt gap (^_e) and C„„{i) (^nn) as a func- 
tion of V/t inside the HI phase, (inset) Haldane-like gap for 
number- conserving excitations after extrapolation to L — > cxd; 
(c) Partition dependence of A2 — A3 compared to Cnn{i) for 
L ^ N ^ 100, for two MI cases. 



A{La) = A2 — A3 in the MI phase. Our analysis shows 
that A{La) also decays exponentially, with a characteris- 
tic length identical to that characterizing the exponential 
decay of C„„(i) (Fig. [UJc)), and hence inversely propor- 
tional to the Mott gap. Note that the relation between 
La dependence of ES degeneracies and density-density 
correlations is also particularly interesting since it may 
allow to evaluate correlation lengths with a minimal ef- 
fort in DMRG simulations. 

By means of the ES we may characterize the phase 
diagram of the EBHM as a function of U/t and V/t. 
We evaluate the largest eigenvalues Aj(i^) of pA for 
different La values, and define Xj = j Sla=i 
We then employ the four largest eigenvalues to define 

= Xf — X2 + — Xj , which characterizes the degen- 
eracy of the entanglement, being zero for perfect double- 
degeneracy of the ES. Proceeding in this way for different 
values of U/t and V/t we obtain the results depicted in 
Fig. [3l Note that the figure shows clearly the appear- 
ance of different regions, although a simultaneous anal- 
ysis of MI and HI demands the use of N = L particles, 
which prevents the destruction or polarization of the edge 
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FIG. 3: Density plot of ( (see text) in the plane U/t, V/t, 
for L = N = 61. In the diagram we identify different phases, 
including MI, HI, CDW, SF and SS. The circles indicate the 
phase boundaries in a 5 = 1 spin chain. 



systems is surprisingly linked to density-density corre- 
lations in both MI and HI, providing as well interest- 
ing insights about other phases. Finally, we note that 
our findings concerning block symmetry and partition- 
dependence are expected to be applicable to other mod- 
els, for both locally- and non-locally-ordered phases. 

We thank T. Vekua and E. Orignac for discussions, and 
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Engineering and Space- Time Research QUEST. 
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FIG. 4: Entanglement spectrum as a function of the partition 
La for L = 100, U/t = 0.5 and V/t = 3.. 



states, and as a result the ES in the HI region is not fully 
degenerated. For U/t > 2 the boundaries are in very 
good agreement with our results obtained from correla- 
tion functions. In Fig. |3]we depicted as well the results 
expected for a S' = 1 chain for the boundaries between 
Haldane, large-D and Neel phases, which are in good 
agreement with the results for C for U /t> 2. For U/t < 2 
it is difficult to get a reliable information about the decay 
of correlations. However, the ES provides interesting in- 
formation about the system even at low U. For small U/t 
a clear separate region (the SF region) occurs for small 
V/t. Moreover, we find that at low U the region of low C 
significantly extends inside the CDW region. We ide ntify 
this region with the apperance of a supersolid phase [27[ . 
Interestingly, this region presents a peculiar La depen- 
dence of the ES (Fig. |4]) characterized by a a periodic 
modulation with a progressively shorter period for grow- 
ing V/t when approaching the CDW phase. We have 
found these intriguing dependence also for SF phases at 
low V/t with non-commesurate fillings. 

In summary, the block structure of the ES and its parti- 
tion dependence in finite systems provide interesting, and 
somewhat unexpected, insights about correlated quan- 
tum systems. This has been illustrated for a ID EBHM, 
for which we have shown that block symmetry provides 
an intuitive understanding about spectral degeneracy in 
the HI, which remarkably persists when low on-site re- 
pulsion breaks particle-hole symmetry. Our results show 
as well that the partition-dependence of the ES in finite 
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